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ABSTRACT 

Star pair latitude observations form the basis for the pole positions 
reported by the Intematicnal Polar Motion Service (IPMS). The IPMS 
processes these observations to produce a mean pole position. However, 
the time series of rav observations contains high frequency information 
which is lost in the calculation of the mean pole. In this study, 2931 
star pair observations are analyzed. A possible large excitation at one 
cycle per solar day is observed. The average power level in the fre- 
quency band of the tesseral tides is seen to be high, although the peaks 
do not occur at the expected tidal frequencies. 
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ANALYSIS OF STAR PAIR LATITUDES 


Polar motion and the variation of latitude have usually been analyzed in the 
imig period domain, primarily for the study of the Chandler wd)ble. However, 
a great volume of data mcists at the high frequencies. The measurement tech- 
nique utilized by the International Polar Motion Service (IPM^ is based on the 
viewing of eightemi star pairs every night by each of five observatories [ see 
Melchior, 1957] . These raw star pair observations are usually n^lected in 
favor of the reported pole positions, although the actual observations are pub- 
lished in the annual reports of the IPMS. 

In this analysis, 61 days of data from 7 May 1964 to 6 July 1964 were 
utilized. These dates were chosen because all five observatories had a large 
number of successful observing sessions. During this period, 2931 star pairs 
were observed by the five stations: Gaithersburg, 691; Ukiah, 724; Mizusawa, 
347; Carloforte, 505; and Kitab, 664. Tl.e data are reasonably distributed among 
the five. 

Prior to studying the power distribution, the data was preprocessed. The 
first step was to remove the average station latitudes. The second step was to 
remove the latitude biases introduced by such errors as star catalc^e declina- 
tion errors. This second step was done by removing the average star pair lati- 
tudes, with the average being taken over all five stations. The biases ranged 
from -0'.'45 to 4-0'.'78 and are listed in Table 1, where the grouo and pair numbers 
are taken from the IPMS system. These biases are analogous to the declination 
corrections, AS, published by the IPMS [ Yumi, 1966, Table 12) , and their values 
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Table 1 


Star pair biases determirod from 61 days oi data. 



for 1964 are also presented in Table 1. The last columu in the table gives the 
difference between the calculated bias and the IPMS decli’^ttiion correction. 
This difference is quite stable at approximately 0'.'14 for all the star pairs. It 
was expected that this difference column would be distributed around O'.'OO, and 
the reason for the variation is not known. 
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Following the pr^roceMlng* the residual latitudes were studied to deter** 
mine udiether the sample had any abnormal error distributi<Nis which might have 
invalidated the analysis. The residual latitudes were found to follow a normal 
distrilRiti<m with a standard deviation of which means that each star pair 
observation has a standard error <d this value. It was concluded that the data 
sample is a reasonable me for statistical analysis. 

The data was now analysed for the power distribution. The time series was 
very irregular, with the spacing from me point to the next varying from 11 min- 
utes to 15 hours. Therdbre, a least squares technique was used. At each fre- 
quency. an amplitude phase were determined for a circularly polarized 
sinusoid. With the frequency steps made small compared to the inverse of the 
data length, the result was a smoothly varying amplitude and phase as a bmction 
of frequency. This type of analysis neglects the correlations which mig^t exist 
between the various physical processes in the variation of latitude. The ampli- 
tude was then squared to generate the power spectrum in Figure 1. 

The typical magnitude for the standard deviation of the amplitude was 0'.'006. 
As was expected, this number is approximately equal to 01'21 (the standard devi- 
ation of the residual analysis above) divided by the square root of the number of 
star pair observations in the data set. However, this O'.'OOS is not a measure of 
the noise level in the data. It is a measure of how well the amplitude is deter- 
mined in the least squares calculation. Since the data time series is not equally 
spaced, statistical tests which might have identified the noise level could not be 
applied. 
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FREQUENCY (Cycles per Sidereal Day) 

Figure 1. Power spectrum generated with 61 days of star pair 
latitude data. 

The largest feature in Figure 1 is the doublet at -2.00 cycles per sidereal 
day (cpd). This feature is not really a doublet. The optical observing technique 
based on meridian passage of a star allows the star catalogue biases to create 
a "blind spot" in the spectrum at integral multiples of 1.00 cpd. This large 
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feature at -2.00 cpd is most likely a single broad peak split by the blind spot, 
m faot, the structure at -1.00 and -fl.OO cpd can be generated by assuming a 
broad peak centered at one cycle per solar day and split by a blind spot at one 
cycle per sidereal day. This is shown in Figure 2 where the line represents 
the actual power spectrum and the crosses represent the theoretical model. 


The theoretical model takes a damped harmonic oscillator line shape 

2 

P X 

F(o,) = 






centered on the solar frequency 0.997 or -0.997 cycles per sidereal day and 
multiplies it by a simple V-shaped filter envelope centered on 1.000 cycles per 
sidereal day 


H(co) 


1. - O) 

0.012 


0.988 < c; < 1.012 


- 1. otherwise 

The best fit indicated by the crosses in Figure 2 is achieved by the following 
parameter values 

Negative Peak P = 26. x O'.'Ol^ 

Q = 130. 

Positive Peak P = 15. x O'.'Ol^ 

Q = 40. 
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POWER (o'.'or-) 



The positive peak In Figure 2 appears 1 1 be shifted to the left according to 
the frequency scale. However, a better local frequency marker Is the minimum 
which Is known to occur at precisely ■•■1.000 cpd. This minimum is also shifted. 
Using the minimum to calibrate the frequency scale causes much of the apparent 
frequency shift to dls{q>pear. 

The general envelope of the spectrum In Figure 1 indicates an interesting 
feature. In the region from -1.2 to -0.8 cpd the spectrum envelope is higher 
than it is in most of the other regions. This frequency band is the location of 
the tesseral constituents of the solid earth rnd ocean tides . These tesseral 
constituents are the only ones in which the perturbing tidal potential has the 
geometry necessary to cause motion in the pole position. Thus, if any tidally 
generate polar motion were present, it would be in this frequency band. How- 
ever, the peaks in this band are not located at the tidal frequencies. 

The primary value of Figure 1 is that it highlights the deficiencies of pre- 
vious high frequency analyses of the latitude variation. These previous analyses 
(for example, Spencer-Jones ( x939J , Markowitz and Bestul [1941] ) sought to 
determine certain elastic parameters of the earth by 'determining the amplitude 
and phase of the latitude variation at precisely the tidal frequencies. Implicit in 
this type of study is the assumption that a local maximum occurs at the tidal 
frequency and that this peak is associated with the tidal potential rather than 
another physical phenomenon. 

In Figure 1, the peaks at the tidal frequencies are conspicuous by their 
absence. If this analysis had proceeded in the previous manner, the results for 
the Oj and Mj tides would have been 0','014 ± 0'.'006 and O'.'OOl ± 01'006, in the 
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retrograde amse respeottvely> we|(1 01*003 ± O'.'OOO and 0V008 * 01*006, in the 
positive sense respectively. Noting that the observing stations are all at the 
same latitude, these values could then have been associated with tidal deflec- 
tions of the vertical and tidally Induced shifts in the principal axes of Che earth. 

When these frequencies are seen in Figure 1 it becomes apparent that this 
simplistic view is pot adequate. The major peaks are not at the tidal frequencies. 
There are a nuxnber of peaks in this region with amplitudes larger than 3. x 
O'.’Ol^. The periodic phenomena which generate these peaks are not understood. 
The only possible tidal peak is the small one at the Oi frequency. When viewed 
with the nearby peaks it becomes impossible to assert that this peak is due to 
the tidal potential rather than some phenomenon associated with the larger 
peaks. Until the major features of the spectrum of the latitude variation are 
understood, conclusions based on the minor features sh6uld be viewed with 
suspicion. 

Hi^ frequency analysis of the variarion of latitude opens up for study a new 
region of the dynamical processes of the earth. The relatively new techniques 
of satellite laser ranging and long baseline radio interferometry offer the oppor- 
tunity for an independent data source and subsequent comparison with Figure 1. 
Present observing schedules are not geared toward providing the high frequency 
data and would have to be altered to at least three latitude measurements per 
day. These techniques are not based on the meridian passage of a star and there- 
fore do not have the blind spot at one cycle per sidereal day. And if the data can 
be acquired as an uiUiroken series of equally spaced data points, then more 
sophisticated statistical tests would be able to be applied. 
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